
* OPEN STATA OUTPUT FILE LOG *

log using "C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.APPENDIX 3.ONLY LONG PANELS.06-02-2022.smcl" 
   
   

use "C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.Post-Estimation.06-02-2022.dta", replace 




xtset a_id year, yearly

  

**** 2022 PAR Data Replication [06/02/2022] "UNDERSTANDING ORGANIZATIONAL SUSCEPTIBILITY TO AGE DISCIMINRATION WITHIN THE U.S. FEDERAL WORKFORCE" [KRAUSE & PARK] ****



 
**** APPENDIX SECTION 3: PANEL SENSITIVITY ANALYSES -- INCLUSION OF ONLY AGENCIES WITH AT LEAST SEVEN ANNUAL OBSERVATIONS [T-BAR = 6.9]  **** 

   
   
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
   
*** GENERATE IDENTIFER FOR PANELS OF VARYING SIZES: FULL PANELS [T=10] = 2; WEAKLY UNBALANCED PANELS [9 =< T >= 7: ABOVE THE MEAN ANUAL OBSERVATIONS = 6.9]; SHORTER PANELS [T < 7] ***

*** 38 agencies contain full panels of observations [T= 10] ***
gen suffpaneldum = 2 if a_id==1 | a_id==8 | a_id==9 | a_id==10 |  a_id== 35 | a_id==38 | a_id==40 | a_id==44 |  a_id== 50 | a_id==51 | a_id==53 | a_id==54 | a_id==55 | a_id==56 | a_id==57 | a_id==58 |  a_id== 59 | a_id==61 | a_id==62 | a_id==65 |  a_id==66 | a_id==75 | a_id==76 | a_id==77 | a_id==79 | a_id==80 | a_id==81 | a_id==83 | a_id==91 | a_id==92 | a_id==96 | a_id==102 |  a_id==112 | a_id==115 | a_id==123 | a_id==124 |  a_id== 127 | a_id==130

*** 43 agencies contain 7-9 observations per panel [9 <= T >= 7] ***

replace suffpaneldum = 1 if a_id==6 | a_id==7 | a_id==12 | a_id==13 | a_id== 14 | a_id==31 | a_id==32 | a_id==33 |  a_id==34 | a_id==37 | a_id==39 | a_id==41 | a_id==42 | a_id==43 | a_id==48 | a_id==49 | a_id==52 | a_id==67 | a_id==68 | a_id==69 | a_id==70 | a_id==71 | a_id==72 | a_id==73 | a_id==74 | a_id==82 | a_id==85 | a_id==87 | a_id==94 | a_id==95 | a_id==97 | a_id==100 | a_id==101 | a_id==103 | a_id==105 | a_id==108 | a_id== 111 | a_id==116 | a_id==117 | a_id==120 | a_id==121 | a_id==125 | a_id==129

*** 49 agencies contain Less than 7 observations per panel [T < 7] ***
replace suffpaneldum = 0 if suffpaneldum==.

tab suffpaneldum 







*** TABLE A3.2: REGRESSION MODEL TABLE PREDICTING VARIATIONS IN AGE DISCRIMINATION FORMAL COMPLAINTS IN U.S. FEDERAL AGENCIES ***





** MODEL A3.1: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- REDUCED MODEL: ONLY CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if suffpaneldum!=0 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio_40over_suplb  ratio_40over_nonsuplb ratio_55over_totallb) 
*
estat ic


*
*
*
*
*

** MODEL A3.2: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem   nonprof40over_tr_lb   politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb  if suffpaneldum!=0, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio_40over_suplb  ratio_40over_nonsuplb    orgjustice_sem  nonprof40over_tr_lb   politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic







** MODEL A3.3: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] REDUCED MODEL: ONLY CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio40suplb_nonsuplb   lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if suffpaneldum!=0 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio40suplb_nonsuplb ratio_55over_totallb) 
*
estat ic

*
*
*
*


** MODEL A3.4: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio40suplb_nonsuplb    orgjustice_sem    nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if suffpaneldum!=0 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio40suplb_nonsuplb   orgjustice_sem  nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic




log close
